Two-band superconductors: Hidden criticality deep in the superconducting state 
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We sliow that two-band superconductors harbor hidden criticality deep in the superconducting 
state, stemming from the critical temperature of the weaker band taken as an independent sys- 
tem. For sufficiently small interband coupling 7 the coherence length of the weaker band exhibits 
a remarkable deviation from the conventional monotonic increase with temperature, namely, a pro- 
nounced peak close to the hidden critical point. The magnitude of the peak scales as oc 7"^, with 
the Landau critical exponent ^ = \, the same as found for the mean-field critical behavior with 
respect to the source field in ferromagnets and ferroelectrics. Here reported hidden criticality of 
multi-band superconductors can be experimentally observed by, e.g., imaging of the variations of 
the vortex core in a broader temperature range. Similar effects are expected for the superconducting 
multilayers. 

PACS numbers: 74.40.Kb, 05.70.Jk, 74.70.Xa, 74.20.De 



Critical phenomena |l| constitute one of the most im- 
portant aspects of the physics of complex systems. The 
classical scenario of both thermal and quantum criticality 
involves the ordered phase induced by a breakdown of a 
basic symmetry and the disordered phase appearing due 
to the restoration of this symmetry. These two phases are 
separated by a critical point, i.e., a second-order phase 
transition accompanied by critical phenomena, e.g., the 
Curie point, the superconducting-to-normal state transi- 
tion, the metal-insulator transition etc. (see, e.g., [lll2|). 
However, possible realizations of the critical behavior are 
not restricted to this standard picture. The present Let- 
ter reports a fascinating example when in addition to the 
standard critical behavior, the system is affected by hid- 
den criticality deep in the ordered phase. This is the case 
of a two-band superconductor [3| (for recent activity see, 
e.g., y42l) where the interband coupling measures the 
proximity to the hidden critical point and plays the role 
of the source field governing the hidden criticality. 

The family of multi-band superconducting materials is 
characterized by the presence of multiple sheets of the 
Fermi surface (bands). In this case the superconducting 
properties are controlled by a set of different band con- 
densates [3[ . Due to the presence of a nonzero interband 
coupling, such condensates are not independent. As a 
result, a band order parameter, i.e., the measure of the 
condensate in a given band, is not simply proportional 
to the Cooper-pair amplitude in this band but involves a 
sum over all band pairing amplitudes, each multiplied by 
a specific coefficient, i.e., the coupling constant gij = gji. 
In a two-band superconductor, when pairing of electrons 
between bands is negligible (e.g., due to symmetry rea- 
sons), the order parameter Ai(x) reads 
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lous average of the field operators with the spin projec- 
tion up and down for the singlet s-wave pairing) . The sys- 
tem described by Eq. ([T]) has a unique critical tempera- 
ture Tc for any nonzero interband coupling 912. However, 
when gi2 = 0, one deals theoretically with two indepen- 
dent superconducting condensates with two different crit- 
ical temperatures Td (for the stronger band) and Tc2 (for 
the weaker band), and the basic symmetry changes from 
U(l) to U(l) X U(l). Here the question arises whether 
or not the behavior of the more realistic, weakly coupled 
system is affected by the proximity of decoupled bands. 
Analytical results. The initial step in our study is to 
get analytical information about a two-gap system with 
weakly coupled components. With this in mind, we ex- 
amine Eq. (m at T = Tc2 by developing a Ginzburg- 
Landau (GL) type of approach for gi2 — > 0. Actu- 
ally, instead of 512 it is more convenient to deal with 
the dimensionless parameter 7 = Ai2/(AiiA22), where 
Xij ~ gijN{Q), and A^(0) is the total density of states 
at the Fermi energy, i.e., iV(0) = Ei^i(O)- The GL 
approach invokes an expansion in powers of the order 
parameter and its spatial derivatives. Such an expan- 
sion for the Cooper-pair amplitude (anomalous Green's 
function) in a two-band superconductor reads 
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with the gauge-invariant gradient D = V — ^ A and the 
parameters (in the clean limit) 
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where i,j enumerate the two bands, and ('0j-|-(x)^j4,(x)) 
is the band-dependent Cooper pair amplitude (anoma- 



Here huJc is the cut-off energy, ^(x) is the Riemann 
zeta-function, F ~ 0.577 is the Euler constant, n^ = 
Ni{0)/N{0), and Vi is the band-dependent Fermi veloc- 
ity. Note however that at T = Tc2 the expansion given by 
Eq. ([2]) holds only for the weaker band in the limit 7^0, 
while it is inappropriate for the stronger band where the 



order parameter does not vanish in this limit. Never- 
theless, Eq. ([2]) provides plenty of information about the 
weaker band. Using Eq. ([T]) and invoking the expansion 
of Eq. © for the weaker band, we find [lO| the following 
7GL equation for A2 at T = Tc2'- 
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with Ai^^^o being the order parameter of the strong 
band in the limit 7 — > 0, and j32 and /C2 given by Eq. ^ 
at T = Tc2- We stress that Eq. Q is exact for A2 in the 
leading order in 7 at T = Tf.2- Note that due to T = Tc,2 
there is no linear term in A2 in Eq. (j4|), in contrast to 
the ordinary GL theory. Despite the absence of the linear 
term in A2 , a stable solution to Eq. ^ exists due to the 
presence of the source term 7Ai_^_yo, which reflects the 
Josephson-like coupling between the bands. This solution 
can be obtained only after solving the proper formalism 
for the stronger band at 7 = 0. Nevertheless, qualitative 
information about A2 as a function of 7 can be found 
from a general analysis of Eq. ^ . In the simplest case 
of a spatially uniform system we obtain for 7 — >■ and 



A2 = ('2^i£^y^' ^A2 _ 7''/' / Ai,^^o y/3 
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Based on Eq. ([S]) and assuming a similar general depen- 
dence of A2 on 7, we can analyze the solution of Eq. ([4]) 
in the presence of a spatially nonuniform condensate. We 
find for 7 -^ and T = Tc2 



A2 0C7^/^ ^oc7-^/^ 6oc7-^/^ 

07 



(6) 



where ^i denotes the band-dependent healing (coherence) 
length and the asymptotic behavior for ^2 is found from 
V^A2 ex 5^^A2 ex 7. Equations ([5|) and ^ therefore 
reveal the critical behavior (see, e.g., Ref. Ill) of the order 
parameter in the weaker band oc 7^' , the corresponding 
susceptibility ex ^^/°~^ and the healing length (propor- 
tional to the correlation radius of fluctuations Tc) oc 7"^, 
with the Landau mean-field critical exponents (5 = 3 and 
/i = i. As seen, the interband coupling 7 can be in- 
terpreted as a source field governing the hidden critical- 
ity in two- band superconductors [11|. This is analogous 
to the source fields in criticality of e.g. ferromagncts 
(magnetic field) or ferroelectrics (electric field), where 
the same Landau exponents are found within the mean 
field. Here we note that A12 is not easily tunable in two- 
band superconductors, which is different as compared to 
the source field in ferromagnets and ferroelectrics. How- 
ever, the coexistence of two weakly coupled order param- 
eters was also achieved in layered mesoscopic rings [121, 
where the interlayer coupling can be varied by chang- 
ing the distance between two participating layers. The 
physics of the coherent phenomena in such a two-layer 



system is essentially the same as in two-band supercon- 
ductors J12lll3|. 

Equation ([2]) enables us to further examine the de- 
pendence of A2 on 7 for T > Tc2- When 7 — )► and 
Tc2<T < Tci we find (for bulk, see 0) 



A2 = 7Ai^-,_>o ?^2 ln(T/Tc2) 



(7) 



where Td is the critical temperature of the stronger band 
for 7 = 0. The behavior of A2 changes dramatically as 
compared to Eq. ((S]). The hidden critical point manifests 



here in the fact that the factor 
when T -^ Tc2 
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diverges 



We should note the existence of 
another (analogous) hidden criticality at Td , which is the 
reason why Eq. ([7]) is not valid for T = Td- However, for 
7 — >■ this point is in close vicinity to the overall critical 
temperature T^. and, as a result, it is always strongly 
overshadowed by the usual critical behavior close to Tc- 

Numerical results. In what follows, we perform a full 
numerical study within a microscopic formalism, in a 
broader range of temperatures and interband couplings. 
In addition, we need to confirm the above analytical re- 
sults since Eq. ^ assumes the validity of the gradient 
expansion - which is only correct when the correspond- 
ing healing length goes to infinity when 7 — >■ 0. In other 
words, the important limitation of our 7GL-analysis is 
that we first assumed ^2 — ^ 00 (7 — >■ 0) in order to sub- 
sequently obtain the diverging behavior of ^2 with the 
critical exponent jjl, which may be misleading. 

As an appropriate theoretical formalism, we choose the 
Bogoliubov-dc Gennes (BdG) equations [la] . We empha- 
size that the standard two-band GL formalism cannot be 
used for our purpose, being well justified only near T,.. 
In the present case, the BdG equations read 
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where Ui,y(x), Wi,y(x) and Eiy are the particle-like and 
hole-like wave-functions and the quasiparticle energy, 
with the subscript i enumerating the band and v the set 
of relevant quantum numbers. The single-electron energy 

reads Ti(x) = ^Sfr'^'^ ~ '^'' '^ith "li the band mass (set 
to the free electron mass) and Hi the chemical potential 
measured from the lower edge of the corresponding band. 
Formally, the BdG equations for different band order pa- 
rameters look decoupled but, in fact, they are connected 
through Eq. ([T|) taken with 

(^,t(x)^,^(x)) = Y. ^-(^)C W [1 - 2/(^..)] , (9) 

V 

where j(Eiv) is the Fermi distribution of quasiparticles. 
Inserting Eq. ((9|) into Eq. ([T]), one needs to remedy the 
ultraviolet divergence: though the product Uiy(x)u*^,(x) 
vanishes for large energies, its decay is not sufficiently 
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FIG. 1: (Color online) (a) The weaker-band order parameter A2 versus the temperature T and the interband coupling 7 as 
calculated from the BdG equations for bulk, (b) A2 as function of 7 for temperatures T — 2.64, 2.65, 2.7 and 2.8 K from 
the BdG equations (solid curves); the dashed curve shows the results of Eq. ([5]) for A2 at T ~ Tc2', the dotted curve is the 
dependence given by Eq. at T = 2.7 K. (c) The 7-susceptibility ^2. as a function of T and 7, that diverges at the critical 
point (T,7) = (r,2,0). 



fast to provide convergence of the sum over the relevant 
quantum numbers given by Eq. ^. Following the stan- 
dard cut-off procedure as implemented in other papers 
on two-band superconductors [3|, la, |7|), the cut-off en- 
ergy hwc is used. 

We first investigate a numerical solution of Eqs. ([5]) 
and (O for a bulk superconductor for small 7. Figure [1] 
shows calculated results for the weaker band in a bulk 
superconductor, for the set of parameters given in the 
figure. Note that this particular choice of the parameters 
is not decisive for our conclusions - similar results are 
found for ^i up to ~ 10 eV, typical of metals. Our choice 
oi Hi ~ 10 meV is justified by recent angle-resolved pho- 
toemission experiments (ARPES) on FeSe2;Tei_a:, which 
revealed the presence of multiple bands with extremely 
small Fermi energies (~ 5meV) |16| . 

Figure [T]Ja) shows that the increase in A2 with 7 be- 
comes faster as T approaches Tc2 ~ 2.645 K, which sup- 
ports our expectations. A2 versus 7 in the uniform sys- 
tem is plotted in Fig. [TJb) for selected temperatures, 
where a comparison is made with our analytical results, 
Eqs. (O and ^. As seen, for T > Tc2 and 7 -> the 
weaker-band order parameter A2 is indeed a linear func- 
tion of 7. However, when T approaches Tc2, this linear 
dependence is replaced by oc 7^^^, as confirmed by the 
very good agreement between the numerically obtained 
curves for T = 2.64 K and 2.65 K and the analytic curve 
for T « Tc2 in Fig.fHb). At T = Tc2 the critical behavior 
of A2 as function of 7 leads to a diverging 7-susceptibility, 
as shown numerically in Fig.[TJc). Due to the well-known 
interrelation between the susceptibility and the coherence 
length, one expects that ^2 is also a divergent function of 
7 at T = Tc2- However, to study this feature and check 
the scaling ^2 oc 7"^/'^ found in the 7GL model, we need 
to abandon the uniform case and study a spatially vary- 
ing two-band condensate. 

We chose to investigate the single-vortex solution in 



a superconducting cylinder (with radius R = 300 nm) 
by numerically solving Eq. ([5]), and extract the band- 
dependent healing lengths ^^ from the size of the vortex 
core. When numerically studying a single-vortex solu- 
tion in the cylinder, we follow a procedure similar to 
that of Ref . [iTJ. We neglect the screening of the magnetic 
field, which is justified for extreme type-II superconduc- 
tors. The spatial variation of the order parameters in 
two bands A^ (p) are shown around a vortex in Fig. ^a) 
in units of their bulk zcro-tcmperature values A^^buik, for 
7 = 0.0002 and T = 0. We defined the healing lengths by 
^i(p = £.i) = O.SAi^buik, i.e., £,1 = 14nm and ^2 = 74nni 
for this particular case [l8| . Figure ^h) shows the tem- 
perature dependence of ^2 for different 7 as extracted 
from our numerical results. We point out our main find- 
ing: the existence of a clear peak in ^2 at temperatures 
close to Tc2- Its peak value increases with decreasing 7 
and almost ideally follows the scaling •y"^'^ found within 
the 7GL model. For sufficiently small interband cou- 
plings, ^2 approaches its independent-band limit shown 
by the dashed curve in Fig. [^Jb). The peak in £,2{T) also 
marks the temperature range where the most pronounced 
difference is found between the two healing lengths in 
the two bands [shown in Fig. [5Jc) for 7 < 0.2]. When 
approaching Tc, the difference between the spatial pro- 
files of the band-dependent condensates disappears, i.e., 
C2/C1 ~^ 1; in agreement with Rcfs. [7', 'sll. The present 
study shows that the effect of the hidden criticality on 
the vortex core is very pronounced for Tc2 < O.STc- At 
higher temperatures the peak in £2 is overshadowed by 
the usual critical behavior around T^, unless the coupling 
7 is extremely weak. 

As seen in Figs.[2l^b,c), with increasing coupling 7 the 
temperature where ^2 peaks, labeled T*, increases while 
the peak itself decreases (for the chosen parameters, the 
peak is washed out for 7 > 0.2). Regarding the maxi- 
mal peak magnitude, it is limited by the taken sample 
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FIG. 2: (Color online) (a) Ai(p) in the vortex core as calcu- 
lated for T = and 7 = 0.0002 and the definition of the heal- 
ing lengths 5i. (b) The numerically extracted healing length 
of the weaker band ^2 as a function of T/Tc for different 7. 
The dashed curve represents ^2 in the independent-band limit 
7 — >■ 0, and the dotted line marks Tc2 = 2.645 K (for com- 
parison, they are shown versus T/Tc with Tc calculated for 
7 = 0.02). (c) The ratio of the healing lengths ^2/Ci versus 
temperature, (d) The bandgaps versus temperature and the 
definition of the crossover temperature T* (where ^2 peaks) 

extracted from the condition -g^ = 0. 



radius (0.3 /im). The largest value of ^2 in Fig. [J^b) is 
well below this limit and finite size effects do not alter 
our conclusions. It is interesting that T* can be evalu- 
ated without a time-consuming numerical study of the 
vortex solution, being very close to the temperature at 
which 9^A2/9r^ = for a bulk superconductor. This is 
illustrated in Fig. [IJd) for 7 = 0.02, where the inflection 
point of A2(T) at T = 0.33 T^ matches weU with T* = 
0.31 Tc found from the peak position in panel (b). In ei- 
ther case, T* should be seen as the crossover temperature 
between the two regimes: lower temperatures - where A2 
is governed by the properties of the weaker band (for suf- 
ficiently small 7), and higher temperatures - where A2 is 
controlled by the tunneling from the stronger band. 

Finally, we briefly discuss available two-gap materials 
from the point of view of the hidden criticality. Based on 
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TABLE I: The coupling constants, calculated nominal critical 
temperatures Tci and Tc2 of the separate bands, the approx- 
imate crossover temperature T*, the critical temperature Tc 
and the interband coupling 7 for selected two-gap materials. 
The coupling constants are taken from the references in the 
last column of the table. 
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In summary, we demonstrated that the properties of 
a two-band superconductor are affected by two critical 
points. In addition to the ordinary critical temperature 
Tc, there exists a hidden critical point - at the critical 
temperature of the weaker band in the absence of cou- 
pling. Interband coupling 7 controls the proximity to this 
hidden critical point and governs the criticality similarly 
to an external magnetic field for ferromagnetic materi- 
als or an external electric field for ferroelectric systems. 
For weak coupling, the weaker-band healing length ex- 
hibits an atypical temperature dependence, with the well 
pronounced peak close to the hidden critical point. This 
gives rise to a large disparity of healing lengths of the 
two condensates in this temperature region. Such a com- 
petition of length scales in a single material can lead to 
significant new physics, and may be closely related to the 
recently observed exotic behavior of vortex matter [J, |5[ . 
Direct measurement of the weaker-band coherence length 
can be realized experimentally, as was demonstrated for 
■K band in MgBj [26|. We suggest that more experimen- 
tal work should be done on recent multiband materials 
close to their hidden critical point, where evidence for 
criticality can be found through unusual thermal proper- 
ties [27| , resurgence of fluctuations deep in the supercon- 
ducting state or anomalous cusps in superfluid density 
as a function of temperature [23l |24| . Further experi- 
mental advancements are facilitated by direct similari- 
ties of two-band superconductors with superconducting 
bilayers, where the atypical temperature behavior of the 
coherence length can be observed as a function of an 
alterable interlayer coupling. Some signatures of such 
behavior are already visible in Ref. |l2 . 

This work was supported by the Flemish Science Foun- 
dation (FWO-Vl). Useful discussions with A. V. Vagov 
are acknowledged. 



the coupling constants available from Refs. 
estimate the dimensionless coupling strength 7 for vari- 
ous two-gap superconductors in Table H] [251 . According 
to our calculations the best candidates to observe the 
effects associated with the hidden criticality seem to be 
VsSi and FeSci-a, due to their particularly low values of 
coupling 7 and low T* and Tc2 in comparison with the 

bulk Tc- This is supported by the experimentally ob- It is convenient to recast the self-consistency equation 

served anomalies close to T* in these two materials, see [i.e., Eq. (1) of the article] for a two-band superconductor 



SUPPLEMENTARY MATERIAL 



Derivation of the 7GL equation 



in the form 

■ Ai(x) ■ 
A2(x) 



All Ai2\ /i?i[Ai(x)] 



A21=Ai2 A22y I i?2[A2(x)] 



, (10) 



where Ri[Ai{x)] is a functional of Ai(x) related to the 
anomalous Green function as [see Eq. (2) in the article] 



R^[\i^)] 



(V'»t(x)^»4(x)) 
N{0) 
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and Xij = Xji = N{Q)gij, with A^(0) the total density of 
states. Then, Eq. (|10l) can further be rearranged as 



i?i[Ai(x)]\ _ 1 / A22 -Ai2\ /Ai(x 



i?2[A2(x)]; jy 1^-Ai2 All J VA2(X 



, (12) 



with 77 = A11A22 — A^2- Based on Eq. (2) of the article 
and Eq. ([T2|) given above, one finds 



f^i-a2)A2+/32A2|A2r-/C2D2A2 ^Ai =0, (13) 

where a2i P2 and /C2 are defined in Eq. (3) of the article. 
Now, assuming that 

r = Te2 = ^fta;ce-^/("=^==^ 

with n.i = N.i{Q)/N{0) the partial density of states, and 
keeping only the terms linear in A12, Eq. (J13p yields 



/32A2|A2r-/C2D2A2 



Ai 



AiiA 



22 



-Ai,Ai,^o = 0, (14) 



with (32 and IC2 taken at T = Tc2, which is the 7GL equa- 
tion [7 = Ai2/(AiiA22)] given by Eq. (4) of the article. 



Derivation of Eq. (7) in the article 



Let us again consider Eq. (|T3l) but now for Tc2 < T < 
Tci and in the homogeneous case. In this temperature 
domain A2 — > for A12 — > and Ai^x^^^q ^ 0. From 
Eq. dni) we have 



A12 



All /2e ftWf 

712 In ( ^7- 

r] \ TTl 



A2 + n2^A2|A2|^ = — Ai. 

(15) 

Note that the expansion given by Eq. (2) in the article 
and therefore also Eq. P^ given above is not justified 
for a spatially nonuniform solution at Tc2 < T < Td 
because the coherence length ^2 associated with a spatial 
variation of the condensate in a weaker band does not 
diverge as A12 goes to zero. In other words, one cannot 
invoke the gradient expansion in the Gor'kov derivation 
of Eq. (2) in the article. However, for homogeneous case 
Eq. (2) of the article is simply an expansion in powers of 



the order parameter A2 and, so, is fully correct. Keeping 
only terms linear in A12, from Eq. (|15|) one finds 



772 ln(T/T,2)A2 



A12 



AllA; 



-Ai, 



22 



(16) 



which is Eq. (7) in the article. 

We remark that at first sight Eq. ([TB]) prescribes that 
A2 is infinite when T — >■ Tc2 and A12 = const. This is 
not true because the validity of the expansion in powers 
of A12 given in Eq. ([T51) requires that A12 < ^{^"^'{T), 
and X12 (T) -^ for T -^ Tc2- Hence, A12 can not be 
considered as constant when using Eq. ([T6| for T — !• Tc2- 
In particular, during the derivation of Eq. (|16p we first 
assumed that 



All /2e hojr 

722 hi 



ttT 
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c2 



^12 
A11A22 



+ OiXt2) (17) 



and, then, we ignored contributions of the order A^2 ^^^ 
higher. However, when the first and second terms in 
the right-hand-side of Eq. ([T7)) become of the same order 
of magnitude, then the contributions oc A^2 can not be 
neglected any more. This makes it possible to estimate 



A 



(lim) 
12 



(T) as 



x\T'iT)^\22 Aii7i2in(r/r,2) 



1/2 



(18) 



which proves that Xf2^\T) -^ when T ^ Tc2. No- 
tice that Eq. ^W^ does not hold for T > T^, where 
Ai >,j2_j.o = 0. Formally, it gives A2 = which means 
that A2 is not linear in A12 for small interband cou- 
plings any longer. This is directly related to the fact that 
T = Tci is one more hidden critical point for the system 
of interest, however it is now associated with a stronger 
band. As pointed out in the article, the hidden criticality 
associated with Td is overshadowed by the usual critical 
behavior around Tr. 
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